SPH based numerical treatment of the interfacial interaction of flow with porous media

In this paper, the macroscopic equations of mass and momentum are developed and discretised based on the Smoothed Particle Hydrodynamics (SPH) formulation for the interaction at an interface of flow with porous media. The theoretical background of flow through porous media is investigated in order to highlight the key constraints which should be satisfied, particularly at the interface between the porous media flow and the overlying free flow. The study aims to investigate the derivation of the porous flow equations, computation of the porosity, and treatment of the interfacial boundary layer. It addresses weak assumptions that are commonly adopted for interfacial flow simulation in particle‐based methods. As support to the theoretical analysis, a 2D weakly compressible SPH (WCSPH) model is developed based on the proposed interfacial treatment. The equations in this model are written in terms of the intrinsic averages and in the Lagrangian form. The effect of particle volume change due to the spatial change of porosity is taken into account and the extra stress terms in the momentum equation are approximated by using Ergun's equation and the Sub‐Particle Scale (SPS) model to represent the drag and turbulence effects, respectively. Four benchmark test cases covering a range of flow scenarios are simulated to examine the influence of the porous boundary on the internal, interface and external flow. The capacity of the modified SPH model to predict velocity distributions and water surface behaviour is fully examined with a focus on the flow conditions at the interfacial boundary between the overlying free flow and the underlying porous media

A parameter-free total Lagrangian smooth particle hydrodynamics algorithm applied to problems with free surfaces

This paper presents a new Smooth Particle Hydrodynamics computational framework for the solution of inviscid free surface flow problems. The formulation is based on the Total Lagrangian description of a system of first-order conservation laws written in terms of the linear momentum and the Jacobian of the deformation. One of the aims of this paper is to explore the use of Total Lagrangian description in the case of large deformations but without topological changes. In this case, the evaluation of spatial integrals is carried out with respect to the initial undeformed configuration, yielding an extremely efficient formulation where the need for continuous particle neighbouring search is completely circumvented. To guarantee stability from the SPH discretisation point of view, consistently derived Riemann-based numerical dissipation is suitably introduced where global numerical entropy production is demonstrated via a novel technique in terms of the time rate of the Hamiltonian of the system. Since the kernel derivatives presented in this work are fixed in the reference configuration, the non-physical clumping mechanism is completely removed. To fulfil conservation of the global angular momentum, a posteriori (least-squares) projection procedure is introduced. Finally, a wide spectrum of dedicated prototype problems is thoroughly examined. Through these tests, the SPH methodology overcomes by construction a number of persistent numerical drawbacks (e.g. hour-glassing, pressure instability, global conservation and/or completeness issues) commonly found in SPH literature, without resorting to the use of any ad-hoc user-defined artificial stabilisation parameters. Crucially, the overall SPH algorithm yields equal second order of convergence for both velocities and pressure

Zoning of mucosal phenotype, dysplasia, and telomerase activity measured by telomerase repeat assay protocol in Barrett's esophagus

Glandular dysplasia in Barrett's esophagus may regress spontaneously but can also progress to cancer. The human telomerase RNA template and the human telomerase reverse transcriptase enzyme which do not, of themselves, correlate strongly with telomerase activity, are too often overexpressed in Barrett's dysplasia to predict individual cancer risk. This study relates telomerase activity, mucosal phenotype, and dysplasia in Barrett's esophagus. Biopsies (n = 256) from squamous esophagus, columnar-lined esophagus every 2 cm, esophago-gastric junction, gastric body, and antrum from 32 patients with long-segment Barrett's esophagus were evaluated by telomerase repeat assay protocol (TRAP). Three biopsies for histology (n = 794) were simultaneously taken at each anatomical level. These and all prior and subsequent biopsies (n = 1917) were reviewed for mucosal phenotypes and dysplasia severity. Intestinal-type Barrett's mucosa was present at all levels in Barrett's esophagus. At least one Barrett's biopsy was TRAP(+) in 22 of 32 patients. TRAP positivity of intestinal-type Barrett's mucosa increased distally, possibly as a consequence of mucosal exposure to acid or bile reflux. Native gastric mucosa was rarely TRAP(+) (1/31 corpus, 2/32 antrum), whereas native squamous mucosa usually was TRAP(+) (31/32). Dysplasia almost always involved intestinal-type Barrett's mucosa (85/87; P</p

Numerical study of nonlinear heat transfer from a wavy surface to a high permeability medium with pseudo-spectral and smoothed particle methods

Motivated by petro-chemical geological systems, we consider the natural convection boundary layer flow from a vertical isothermal wavy surface adjacent to a saturated non-Darcian high permeability porous medium. High permeability is considered to represent geologically sparsely packed porous media. Both Darcian drag and Forchheimer inertial drag terms are included in the velocity boundary layer equation. A high permeability medium is considered. We employ a sinusoidal relation for the wavy surface. Using a set of transformations, the momentum and heat conservation equations are converted from an (x, y) coordinate system to an (x,η) dimensionless system. The two-point boundary value problem is then solved numerically with a pseudo-spectral method based on combining the Bellman–Kalaba quasi linearization method with the Chebyschev spectral collocation technique (SQLM). The SQLM computations are demonstrated to achieve excellent correlation with smoothed particle hydrodynamic (SPH) Lagrangian solutions. We study the effect of Darcy number (Da), Forchheimer number (Fs), amplitude wavelength (A) and Prandtl number (Pr) on the velocity and temperature distributions in the regime. Local Nusselt number is also computed for selected cases. The study finds important applications in petroleum engineering and also energy systems exploiting porous media and undulating (wavy) surface geometry. The SQLM algorithm is shown to be exceptionally robust and achieves fast convergence and excellent accuracy in nonlinear heat transfer simulations

SPH study of the evolution of water–water interfaces in dam break flows

The mixing process of upstream and downstream waters in the dam break flow could generate significant ecological impact on the downstream reaches and influence the environmental damages caused by the dam break flood. This is not easily investigated with the analytical and numerical models based on the grid method due to the large deformation of free surface and the water-water interface. In this paper, a weakly compressible Smoothed Particle Hydrodynamics (WCSPH) solver is used to study the advection and mixing process of the water bodies in two-dimensional dam-break flows over a wet bed. The numerical results of the mixing dynamics immediately after the release of the dam water are found to agree satisfactorily with the published experimental and numerical results. Then further investigations are carried out to study the interface development at the later stage of dambreak flows in a long channel. The analyses concentrate on the evolution of the interface at different ratios between the upstream and downstream water depths. The potential capabilities of the mesh-free SPH modelling approach for predicting the detailed development of the water-water interfaces are fully demonstrated.The first author acknowledges the Jafar Studentship during her PhD study at the University of Cambridge. The other authors acknowledge the support of the Major State Basic Research Development Program (973) of China (No. 2013CB036402), Open Fund of the State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University (SKHL1404; SKHL1409), Start-up Grant for the Young Teachers of Sichuan University (2014SCU11056) and National Science and Technology Support Plan (2012BAB0513B0).This is the accepted manuscript. The final version is available at http://link.springer.com/article/10.1007%2Fs11069-015-1726-6

A rarefaction-tracking method for hyperbolic conservation laws

We present a numerical method for scalar conservation laws in one space dimension. The solution is approximated by local similarity solutions. While many commonly used approaches are based on shocks, the presented method uses rarefaction and compression waves. The solution is represented by particles that carry function values and move according to the method of characteristics. Between two neighboring particles, an interpolation is defined by an analytical similarity solution of the conservation law. An interaction of particles represents a collision of characteristics. The resulting shock is resolved by merging particles so that the total area under the function is conserved. The method is variation diminishing, nevertheless, it has no numerical dissipation away from shocks. Although shocks are not explicitly tracked, they can be located accurately. We present numerical examples, and outline specific applications and extensions of the approach.Comment: 21 pages, 7 figures. Similarity 2008 conference proceeding

Associations between maternal thyroid function in pregnancy and child neurodevelopmental outcomes at 20 months in the Seychelles Child Development Study, Nutrition Cohort 2 (SCDS NC2)

Maternal thyroid hormones facilitate optimal foetal neurodevelopment; however, the exact role of the thyroid hormones on specific cognitive outcomes is unknown. The present study aimed to investigate associations between maternal thyroid function and neurodevelopmental outcomes in the Seychelles Child Development Study (SCDS) Nutrition 2 cohort (n 1328). Maternal free thyroid hormones (fT3, fT4 and fTSH) were assessed at 28 weeks’ gestation with a range of child cognitive outcomes analysed at 20 months. Dietary iodine intake was analysed for a subset of women through a Food Frequency Questionnaire. Linear regression analysis was used to test associations between serum concentrations of maternal thyroid hormones and child neurodevelopment outcomes. Thyroid hormones were analysed as continuous data and categorised as quintiles. 95% of mothers had optimal thyroid function based on fTSH concentrations. Overall, the present study shows that maternal thyroid function is not associated with neurodevelopmental outcomes in this high fish-eating population. However, a positive association, using quintiles for fT3, was reported for the Mental Developmental Index, between Q3 v. Q4 (β 0⋅073; P 0⋅043) and for Q3 v. Q5 (β value 0⋅086; P 0⋅018). To conclude, mothers in our cohort, who largely have optimal thyroid function and iodine intakes, appear able to regulate thyroid function throughout pregnancy to meet neurodevelopmental needs. However, it is possible that minor imbalances of fT3, as indicated from our secondary analysis, may impact offspring neurodevelopment. Further investigation of the relationship between maternal thyroid function and infant neurodevelopment is warranted, particularly in populations with different dietary patterns and thereby iodine intakes

Potential application of mesh-free SPH method in turbulent river flows

A comprehensive review has been completed on the simulation of turbulent flow over rough beds using mesh-free particle models. Based on the outcomes of this review, an improved Smoothed Particle Hydrodynamics (SPH) method has been developed for open channel flows over a rough bed, in which a mixing length model is used for modeling the 2D turbulence and a drag force equation is proposed for treating the boundary shear. The proposed model was applied to simulate a depth-limited open channel flow over a rough bed surface. The results of the velocity profile and shear stress distribution show a good agreement with the experimental data and existing analytical solutions. This work reveals that in order to correctly model turbulent open channel flow over a rough bed, the treatment of both flow turbulence and bed roughness effect is equally important

Veterinary dairy herd fertility service provision in seasonal and non-seasonal dairy industries - a comparison

The decline in dairy herd fertility internationally has highlighted the limited impact of traditional veterinary approaches to bovine fertility management. Three questionnaire surveys were conducted at buiatrics conferences attended by veterinary practitioners on veterinary dairy herd fertility services (HFS) in countries with a seasonal (Ireland, 47 respondents) and non-seasonal breeding model (The Netherlands, 44 respondents and Portugal, 31 respondents). Of the 122 respondents, 73 (60%) provided a HFS and 49 (40%) did not. The majority (76%) of all practitioners who responded stated that bovine fertility had declined in their practice clients' herds with inadequate cow management, inadequate nutrition and increased milk yield as the most important putative causes. The type of clients who adopted a herd fertility service were deemed more educated than average (70% of respondents), and/or had fertility problems (58%) and/or large herds (53%). The main components of this service were routine postpartum examinations (95% of respondents), fertility records analysis (75%) and ultrasound pregnancy examinations (69%). The number of planned visits per annum varied between an average of four in Ireland, where breeding is seasonal, and 23 in Portugal, where breeding is year-round. The benefits to both the practitioner and their clients from running a HFS were cited as better fertility, financial rewards and job satisfaction. For practitioners who did not run a HFS the main reasons given were no client demand (55%) and lack of fertility records (33%). Better economic evidence to convince clients of the cost-benefit of such a service was seen as a major constraint to adoption of this service by 67% of practitioners